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Understanding the physical nature of the dark energy which appears to drive the accelerated 
expansion of the unvierse is one of the key problems in physics and cosmology today. This important 
problem is best studied using a variety of mutually complementary approaches. Daly and Djorgovski 
(2003, 2004) proposed a model independent approach to determine a number of important physical 
parameters of the dark energy as functions of redshift directly from the data. Here, we expand this 
method to include the determinations of its potential and kinetic energy as functions of redshift. 
We show that the dark energy potential and kinetic energy may be written as combinations of the 
first and second derivatives of the coordinate distance with respect to redshift. We expand the 
data set to include new supernova measurements, and now use a total of 248 coordinate distances 
that span the redshift range from zero to 1.79. First and second derivatives of the coordinate 
distance are obtained as functions of redshift, and these are combined to determine the potential 
and kinetic energy of the dark energy as functions of redshift. An update on the redshift behavior 
of the dimensionless expansion rate E{z), the acceleration rate q{z), and the dark energy pressure 
p{z), energy density f{z), and equation of state w{z) is also presented. We find that the standard 
Dom = 0.3 and Da = 0.7 model is in an excellent agreement with the data. We also show tentative 
evidence that the Cardassian and Chaplygin gas models in a spatially flat universe do not fit the 
data as well. 


I. INTRODUCTION 

The acceleration of the universe at the present epoch 
can be studied usiM a variety of techniques such as type 
la supernovae [E 0 0) S IS 0. 0 > ^-nd studies of cosmic 
microwave background radiation ( e.g ., I p ) combined with 
large scale structure studies (e.g. |S[i3)) which indicate 
that at the present epoch the universe is expanding at an 
accelerting rate. Studies of powerful FRII radio galaxies 
also indicate that the universe is accelerting [TE [lE US ■ 
To study this in detail, it is important to determine the 
redshift behayior of the expansion and acceleration rates 
of the universe, and the properties and redshift evolution 
of the driver(s) of these rates. 

In this vein, El E3 suggest a model-indepedent 
method of studying the redshift evolution of the expan¬ 
sion and acceleration rates of the universe, and show that 
the acceleration of the universe can be written in terms of 
the first and second derivatives of the coordinate distance 
with respect to redshift. These derivatives are obtained 
directly from the data on coordinate distances using a 
statistically robust numerical technique. The accelera- 
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tion thus derived requires very few assumptions, relying 
only upon the assumptions that the uniyerse is homoge¬ 
neous and isotropic on large scales, and has zero space 
curvature. Indeed, this approach allows a determination 
of the acceleration of the universe that is independent of 
a theory of gravity and of the properties and redshift evo¬ 
lution of the drivers of the expansion and acceleration, 
such as dark energy, dark matter, or other components. 
Complementary model-independent approaches based on 
an integral rather than a differential technique are dis¬ 
cussed by El El El El m iE El m . 

Understanding the properties of the physical driver of 
the acceleration of the universe is of fundamental im¬ 
portance. The driver, commonly referred to as the dark 
energy, and its properties are parameterized in quantities 
such as the potential and kinetic energy, and the pressure, 
energy density, and equation of state. These quantities 
can be expressed as combinations of the first and second 
derivative of the coordinate distance, as discussed by El 

Here, we show that the dark energy potential and ki¬ 
netic energy may be determined as functions of redshift 
by appropriately combining the first and second deriva¬ 
tives of the coordinate distance to sources at different 
redshift. This expansion of the method is described in 
section m The method is applied to an enlarged data 
sample which includes the 71 new Legacy supernovae co¬ 
ordinate distances, presented in section EE and the co- 
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ordinate distances listed in 0 

The work presented here on the dark energy potential, 
V(z), and kinetic, K(z), energy is complementary to the 
work of [3, H and Pa- Saini et al. [3 derive equa¬ 
tions for V(z) and K{z) and use a fitting function for the 
luminosity distance to obtain values and uncertainties for 
the parameters of the htting function; these are then used 
to obtain V{z), K{z), and the equation of state parame¬ 
ter w{z). Simon et al. [3 consider the reconstruction of 
the dark energy potential based on an expansion of the 
dark ener gy p otential in terms of Chebyshev polynomi¬ 
als, while 13 consider the expansion of the quintessence 
potential as a power series of the quintessence field cj). 
13 j and 13 find that the results are consistent with 
those expected if a cosmological constant is driving the 
acceleration of the universe at the present epoch. 


II. THE POTENTIAL AND KINETIC ENERGY 
OF THE DARK ENERGY 

The potential and kinetic energy of the dark energy 
can be determined as functions of redshift using measure¬ 
ments of the coordinate distance. Measurements of the 
coordinate distance obtained using type la supernovae 
and type Ilb radio galaxies allow determinations of the 
first and second derivative of the coordinate distance with 
respect to redshift. For detailed discussion of numerical 
methods used to derived these quantities as well as exten¬ 
sive tests of the methods, see [30- These derivatives 
may be combined to obtain V(z) and K{z), the potential 
and kinetic energy density of the dark energy, which, for 
convenience, are expressed in units of the critical density 
at the current epoch. 

It is well known that V = 0.5(p — P) and K = 0.5(p -I- 
P), where p and P are the energy density and pressure of 
the dark energy. In |I5| we show that both p and P may 
be written in terms of the hrst and second derivatives 
of the coordinate distance (see Eqs. [6] and [7] of [3). 
Combining these, we find that 

f— ) = {vT^ [l + (I+z)(y")(2/')-V3]-0.5Ho^(l+z)3 

\PocJ 

( 1 ) 

and 

f—) =-(I/3)(l+z)(y")(^/')■'-0.5F!o,.(l+^)^ (2) 

\PocJ 

where poc is the critical density at the current epoch, the 
dimensionless coordinate distance y{z) = i7o(aor), Hq is 
Hubble’s constant, (aor) is the coordinate distance to a 
source at redshift z, y' = (dy/dz), and y" = {d^y/dz"^). 

In obtaining Eqs. [[] and [5], it has been assumed 
that: the universe is spatially flat; the kinematics of the 
universe are accurately described by general relativity; 
and two components, the dark energy and non-relativisitc 
matter, are sufficient to account for the kinematics of the 


universe out to redshift of about 2 (see the discussion in 

0 )- 

Thus, the redshift behavior of the potential and kinetic 
energy densities of the dark energy can be constructed 
using the first and second derivatives of the dimension¬ 
less coordinate distance y{z), for the data set described 
below. 


HI. THE DATA SET 

The core of our data set remains the sample of 157 
“Gold” supernovae of 0 and the 20 radio galaxids of 
13 ) both of which are tabulated in . This data set is 
supplemented by 71 new supernvae from the Supernova 
Legacy Survey of 0, which allows the determination of 
dimensionless coordinate distances to 71 additional type 
la supernovae. These are listed in Table I, and were 
obtained using the values and uncertainties of pB listed 
in Table 9 of Q, given the value of Hq = 70 km/s/Mpc 
adopted by 0. Two uncertainties of y are listed: a{y) is 
obtained from the uncertainty in ps listed in Table 9 of 
0 , while the total uncertainty of y (used throughout this 
paper) uriy) reflects the total uncertainty of /is, which is 
obtained by adding in quadrature the uncertainty listed 
in Table 9 to the intrinsic dispersion of 0.13 identified by 
0 . These are listed in Tabled 

The total sample of 248 sources, including 228 super¬ 
novae and 20 radio galaxies, is shown in Fig. ^ We note 
that there are no systematic differences seen among the 
three groups of measurements in the redshift ranges of 
their overlaps. 


IV. RESULTS OBTAINED WITH THE FULL 
SAMPLE OF 248 SOURCES 

Values of y' and y" obtained with the full sample are 
shown in Figs. [ 2 ], and As in 00 , a window func¬ 
tion with Az of 0.6 is used in this paper. Note that the 
first and second derivatives are obtained directly from 
the measured values of y; no assumptions have been 
adopted regarding a theory of gravity, or the form of 
the expansion or acceleration rate of the universe as a 
function of redshift. For comparison, the values of pa¬ 
rameters predicted in a standard lambda cold dark mat¬ 
ter (LCDM) cosmological model, where the kinematics 
of the universe are described by general relativity and 
the primary components of the universe at present are 
a normalized cosmological constant Ha = 0.7 and non- 
relativistic matter Hom = 0.3, are shown on each fig¬ 
ure. In addition, the curve expected in a universe with 
space curvature k is shown in Fig. j2], where values of 
Hom = 0.3, Ha = 0, and H/j = —k/{aoHo)‘^ = 0.7 have 
been assumed. Clearly, this curvature dominated model 
is inconsistent with the results obtained here. 

To compare our results with those expected in a LCDM 
model a mock data set of 248 sources was constructed in 
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TABLE I: Coordinate distances to the 71 Legacy Supernovae 
and their uncertainties. 


name 

z 

y 

^ (y) 

or (y) 

SNLS-03Dlau 

0.504 

0.475 

0.0085 

0.0297 

SNLS-03Dlaw 

0.582 

0.556 

0.0138 

0.0360 

SNLS-03Dlax 

0.496 

0.426 

0.0074 

0.0265 

SNLS-03Dlbp 

0.346 

0.325 

0.0031 

0.0197 

SNLS-03Dlcni 

0.870 

0.822 

0.1140 

0.1242 

SNLS-03D1CO 

0.679 

0.665 

0.0269 

0.0480 

SNLS-03Dlew 

0.868 

0.743 

0.1177 

0.1258 

SNLS-03Dlfc 

0.331 

0.271 

0.0016 

0.0163 

SNLS-03Dlfl 

0.688 

0.562 

0.0127 

0.0360 

SNLS-03Dlfq 

0.800 

0.647 

0.0268 

0.0471 

SNLS-03Dlgt 

0.548 

0.554 

0.0204 

0.0389 

SNLS-03D3af 

0.532 

0.502 

0.0192 

0.0357 

SNLS-03D3aw 

0.449 

0.380 

0.0077 

0.0240 

SNLS-03D3ay 

0.371 

0.338 

0.0047 

0.0207 

SNLS-03D3ba 

0.291 

0.286 

0.0044 

0.0177 

SNLS-03D3bh 

0.249 

0.243 

0.0022 

0.0147 

SNLS-03D3CC 

0.463 

0.417 

0.0065 

0.0258 

SNLS-03D3cd 

0.461 

0.407 

0.0109 

0.0267 

SNLS-03D4ag 

0.285 

0.254 

0.0018 

0.0153 

SNLS-03D4at 

0.633 

0.605 

0.0178 

0.0404 

SNLS-03D4cn 

0.818 

0.653 

0.0914 

0.0994 

SNLS-03D4CX 

0.949 

0.602 

0.0754 

0.0836 

SNLS-03D4cy 

0.927 

0.986 

0.1725 

0.1823 

SNLS-03D4CZ 

0.695 

0.554 

0.0219 

0.0398 

SNLS-03D4dh 

0.627 

0.508 

0.0082 

0.0315 

SNLS-03D4di 

0.905 

0.676 

0.0803 

0.0899 

SNLS-03D4dy 

0.604 

0.463 

0.0062 

0.0284 

SNLS-03D4fd 

0.791 

0.610 

0.0214 

0.0423 

SNLS-03D4gf 

0.581 

0.526 

0.0114 

0.0335 

SNLS-03D4gg 

0.592 

0.477 

0.0198 

0.0347 

SNLS-03D4gl 

0.571 

0.462 

0.0149 

0.0314 

SNLS-04Dlag 

0.557 

0.476 

0.0064 

0.0292 

SNLS-04Dlaj 

0.721 

0.594 

0.0290 

0.0459 

SNLS-04Dlak 

0.526 

0.517 

0.0131 

0.0336 

SNLS-04D2cf 

0.369 

0.338 

0.0025 

0.0204 

SNLS-04D2fp 

0.415 

0.373 

0.0046 

0.0228 

SNLS-04D2fs 

0.357 

0.334 

0.0028 

0.0202 

SNLS-04D2gb 

0.430 

0.370 

0.0065 

0.0231 

SNLS-04D2gc 

0.521 

0.472 

0.0117 

0.0306 

SNLS-04D2gp 

0.707 

0.607 

0.0361 

0.0512 

SNLS-04D2iu 

0.691 

0.587 

0.0368 

0.0509 

SNLS-04D2ja 

0.741 

0.649 

0.0350 

0.0523 

SNLS-04D3CO 

0.620 

0.581 

0.0161 

0.0383 

SNLS-04D3cp 

0.830 

0.973 

0.1556 

0.1661 

SNLS-04D3cy 

0.643 

0.571 

0.0155 

0.0376 

SNLS-04D3dd 

1.010 

0.999 

0.2451 

0.2523 

SNLS-04D3df 

0.470 

0.451 

0.0066 

0.0278 

SNLS-04D3do 

0.610 

0.525 

0.0094 

0.0328 

SNLS-04D3ez 

0.263 

0.253 

0.0015 

0.0152 

SNLS-04D3fk 

0.358 

0.339 

0.0020 

0.0204 

SNLS-04D3fq 

0.730 

0.613 

0.0212 

0.0424 

SNLS-04D3gt 

0.451 

0.411 

0.0057 

0.0253 

SNLS-04D3gx 

0.910 

0.868 

0.1384 

0.1478 

SNLS-04D3hn 

0.552 

0.467 

0.0075 

0.0290 

SNLS-04D3is 

0.710 

0.589 

0.0209 

0.0410 

SNLS-04D3ki 

0.930 

0.930 

0.1841 

0.1924 

SNLS-04D3kr 

0.337 

0.312 

0.0014 

0.0187 

SNLS-04D3ks 

0.752 

0.573 

0.0238 

0.0418 

SNLS-04D31p 

0.983 

0.723 

0.1650 

0.1706 



redshift 


FIG. 1: Coordinate distances to the 71 Legacy and 157 Gold 
supernovae and 20 radio galaxies. 


TABLE II: Table 1 Continued 


name 

z 

y 

(y) 

CTT (y) 

SNLS-04D31U 

0.822 

0.655 

0.0658 

0.0766 

SNLS-04D3nil 

0.950 

0.739 

0.0912 

0.1014 

SNLS-04D3nc 

0.817 

0.690 

0.0807 

0.0907 

SNLS-04D3nh 

0.340 

0.320 

0.0018 

0.0193 

SNLS-04D3nr 

0.960 

0.631 

0.0680 

0.0778 

SNLS-04D3ny 

0.810 

0.706 

0.0978 

0.1065 

SNLS-04D3oe 

0.756 

0.652 

0.0174 

0.0427 

SNLS-04D4an 

0.613 

0.566 

0.0159 

0.0374 

SNLS-04D4bk 

0.840 

0.628 

0.0535 

0.0654 

SNLS-04D4bq 

0.550 

0.473 

0.0122 

0.0308 

SNLS-04D4dni 

0.811 

0.794 

0.0966 

0.1077 

SNLS-04D4dw 

0.961 

0.751 

0.1003 

0.1099 


which the mock sources have the same redshift distribu¬ 
tion and fractional error per point as the empirical data 
but have values of y predicted in a LCDM model with 
Ha = 0.7 and Hom = 0.3. This mock data set was passed 
through the same numerical programs as the true data, 
and the results obtained for y' and y" for both the true 
and mock data sets are shown in Figs. P] and [H]. It is 
quite clear that no bias is introduced by the numerical 
differentiation; the central value of y' and y" output by 
the numerical differentiation lie on those predicted in the 
assumed LCDM model, and the magnitude of the error 
bars at a given redshift have the same magnitude as those 
obtained from the data. 

Values of the dimensionless expansion rate E{z), and 
the deceleration parameter q{z) are shown in Figs. [H] 
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FIG. 2: The first derivative of the coordinate distance with 
respect to redshift. The values for the standard LCDM 
model with Qa ~ 0.7 and Ilom = 0.3 are shown as the 
solid line in this and all subsequent plots. The values for the 
curvature-dominated model described in the text are shown 
with the dotted line. The zero redshift value we measure is 
I/O = 1.025 ± 0.022. The predicted value is 1.000. 


and ini, obtained using equations [5] and [6] of ^3- The 
only assumptions that must be adopted to construct E 
and q from y' and y" are that the universe is homoge¬ 
neous and isotropic on large scales, and is spatially flat 
(see ^3)- Again, the curve expected in a LCDM model 
with riA = 0.7 and Dom = 0.3 is shown in each figure. 
In addition, the curves predicted in two modified gravity 
models is a spatially flat universe are included in Fig. [^. 
These are obtained using the best fit model parameters 
obtained by for the Cardassian model of and 
the generalized Chaplygin gas model of based on the 
model of [^. Clearly, the LCDM model provides a bet¬ 
ter description of the data than do either of the modified 
gravity models. 

Fig. 0 shows the difference between values of E{z) 
and those expected in a LCDM model with Da = 0.7 
and Dom = 0.3. This further illustrates how well Gen¬ 
eral Relativity with a cosmological constant describes the 
data over the very large length scale of greater than about 
10 billion light years. Recall that E{z) is derived from 
the data without having to specify a theory of gravity. 

The results shown for q{z) allow a determination of 
the redshift at which the universe transitions from an ac¬ 
celerating phase to a decelerating phase; we find this to 
be at a redshift of Z't = 0.42 consistent with the 

values quoted by [Tj, [13 and 13 • The upper bound on 



FIG. 3: As in Fig. El but for the mock data set with the same 
redshift and error distribution as the actual data set, with the 
assumed cosmology with Da = 0.7 and Dom = 0.3. The fit 
results for this mock data set are in an excellent agreement 
with the analytical prediction for the assumed cosmology, giv¬ 
ing us a confidence that the numerical method we use is un¬ 
biased and is capable of recovering accurately the underlying 
cosmology. 

this transition redshift is uncertain because of the fluctu¬ 
ations in q{z) which are due to sparse sampling at high 
redshift. In this plot, and in the ones that follow, we do 
not consider these fluctuations at higher redshifts to be 
statistically significant, as they are commensurate with 
our derived 1-cr error bars. 

A comparison with these results and those expected 
in a LCDM model is shown in Figs. [7| and uni, where 
the predicted values of E{z) and q{z) are obtained using 
the mock data set described above. Again, we see that no 
bias has been introduced by the numerical differentiation 
technique, and the data are in very good agreement with 
the standard LCDM model. 

To obtain the pressure, a theory of gravity must be 
specified, and general relativity has been assumed here 
(see eqs. [5] and [6] of M)- Note that if the dark en¬ 
ergy is a cosmological constant, the zero redshift value of 
p = P/ Poc is a measure of Da. The pressure of the dark 
energy, in units of the critical density today, obtained 
with the current data set is shown in Fig. in The zero 
redshift value suggests Da = 0.61 ± 0.08, in excellent 
agreement with the values commonly derived using more 
traditional approaches. 

Values of p{z) expected in a LCDM model are shown 
in Fig. where the predicted values of p{z) are ob¬ 

tained using the mock data set described above. Again, 
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FIG. 4: The second derivative of the coordinate distance with 
respect to redshift. The measured zero redshift value is y'^ = 
—0.55 ±0.10; the value predicted in a standard LCDM model 
is -0.45. 



FIG. 5: As in Fig. H but for the same mock data set. Again, 
the correct assumed cosmology is recovered with a negligible 
bias. 



FIG. 6: The dimensionless expansion parameter E(z). The 
measured zero redshift value is Eq — 0.98 ± 0.02. The pre¬ 
dicted value is 1.00. The solid line shows the prediction of 
the standard Qa = 0.7 and flom = 0.3 model, which is in an 
excellent agreement with the data. The dotted line represents 
predictions of the Gardassian model, and the dashed line rep¬ 
resents the Chaplygin gas model, obtained assuming that the 
universe is spatially flat, both of which seem to fit the data 
systematically less well. 



FIG. 7: Derived values of E{z) for the mock data set. Again, 
the correct cosmology is recovered accurately by the method. 
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redshift 

FIG. 8: Values of E(z) relative to the LCDM model with 
Ha = 0.7 and Qom = 0.3, which has been subtracted from the 
data and the models. This shows how well the LCDM model 
describes E(z) and that the Cardassian model and generalized 
Chaplygin gas model in a spatially flat universe seem to fit 
the data systematically less well. 



redshift 


FIG. 9: The deceleration parameter q{z). The zero redshift 
value is go = —0.46 ±0.08. The predicted value in a standard 
LCDM model is —0.55. Our fits are systematically too high 
by about Icr relative to the standard model. 



FIG. 10: Derived values of the q{z), but for the mock data 
set. Again, the underlying cosmology is recovered, except for 
a slight bias at the high redshift end. 


there is no bias introduced, by the numerical differentia¬ 
tion technique, and the output is in good agreement with 
predictions in a LCDM model. 

The energy density f{z) = pjPoc of the dark energy in 
units of the critical density today can be obtained once 
a value for Dom has been adopted. The value of f{z) ob¬ 
tained assuming Dom = 0.3 is shown in Fig. These 
results are compared with predictions in a LCDM model 
obtained with the mock data set in Fig. m There is 
good agreement between the data and predictions in a 
LCDM model. 

The equation of state of the dark energy, w = p/P 
is shown in Fig. ca Results obtained with the mock 
data set obtained in a LCDM model are shown in Fig. 
nn]. Figs. 11)1 m im [T^ and provide an update on the 
results presented by ( 13 ; the dependence of each quantity 
on the first and second derivatives of the dimensionless 
coordinate distance y with respect to redshift z is given 
in [l5|. 

The values of V and K obtained with the full sample of 
248 sources are shown in Figs. ini and m The results 
are consistent with a dark energy potential energy that is 
constant from a redshift of zero to a redshift of about 0.8, 
as expected if the dark energy is a cosmological constant. 
The zero redshift normalization depends on y'(z = 0) 
and y"{z = 0) (see Eq. [[]). The value of y” at z=0 is 
—0.55 ± 0.1, while that of y' is 1.03 ± 0.02; these cause 
the value of V(z=0) to be low relative to a LCDM model 
with Dom = 0.3. When the number of low-redshift type 
la supernovae is substantially increased, we expect to be 
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redshift 

FIG. 11: The pressure of the dark energy in units of the 
critical density today. The zero redshift value of po = —0.61± 
0.08. 

able to determine whether the normalization of V differs 
from that expected in a standard LCDM model. These 
are compared with predictions in a LCDM model using 
the mock data set, as shown in Figs. m and EOI; a 
LCDM model provides a good description of the data. 

We do not understand completely the small system¬ 
atic difference between our evaluations of some of the 
observed trends from the standard LCDM model: our 
values being systematically too high for the q{z), P{z), 
and w{z), and too low for V(z), by about 1-cr through¬ 
out. We are inclined to interpret this as the inherent 
limitation of the current data set, rather than as a real 
physical effect. 


V. SUMMARY AND CONCLUSIONS 

The redshift behavior of the potential and kinetic en¬ 
ergy of the dark energy can be determined directly from 
observed quantities without having to assume an a priori 
functional form for these quantities. They depend upon 
the first and second derivatives of the coordinate distance 
with respect to redshift, and on the zero redshift value of 
the mean mass density in non-relativistic matter. These 
quantities, y', y", V(z), and K{z) were obtained using a 
sample of 248 sources including 228 type la supernovae 
and 20 type Ilb radio galaxies. In all cases, the results 
are consistent with predictions of a LCDM model out to 
a redshift of about one. 

The data may also be used to determine the decel- 



FIG. 12: Values of the P{z), but for the mock data set. The 
assumed cosmology is recovered accurately. 



redshift 


FIG. 13: The energy density of the dark energy in units of 
the critical density today, f{z). The zero redshift value is 
fo = 0.64 ± 0.04, again in an excellent agreement with the 
more traditional determinations of Da. 





























































































FIG. 14: The values of f{z) for the mock data set. The 
assumed value of 0.7 is recovered accurately. 


eration parameter q(z) and the dimensionless expansion 
parameter E(z) directly from y' and y"] these determi¬ 
nations of q and E assume only that the universe is ho¬ 
mogeneous and isotropic on large scales, and is spatially 
flat, and are independent of any assumptions regarding a 
theory of gravity or the properties of the dark energy, as 
discussed in detail by [ij, [ 13 . We find that the universe 
transitions from an accelerating phase to a decelerating 
phase at a redshift of zt = 0-42±g Qg. The fact that 
and y" match predictions in a LCDM model to about 
one sigma or better to a redshift of about 1.5 means that 
General Relativity with a cosmological constant provides 
an accurate description of the data on scales of about 10 
billion light years. The correct explanation of the dark 
energy must be able to account for the observed values 
of y{z), y'{z) and y"{z), or, equivalently, y{z), E{z), and 
q{z). 

By adopting general relativity as the correct theory 
of gravity, the pressure of the dark energy can be ob¬ 
tained as a function of redshift. As shown by this 
depends only upon y' and y"; it is independent of as¬ 
sumptions regarding the properties of the dark energy 
and independent of Oom. If the universe is dominated by 
a cosmological constant at the present epoch, the zero 
redshift value of p yields a new method of determining 
riA, and the value obtained here is 0.61±0.08. This is re¬ 
markably close to values obtained with more traditional 
approaches. 

The energy density and equation of state of the dark 
energy may be obtained as functions of redshift if the 
value of flom is known. These were determined assuming 



FIG. 15: The dark energy equation of state parameter w{z). 
The zero redshift value is wo = —0.95±0.07, whereas ui = — 1 
is the theoretical value for the cosmological constant model. 



2 _^_I_I_I_I_I_^_I_I_I_I_I_I_^_I_ 1 1 IIINNNIIINNI IIIIIIIIIIIIIIIIIII 

0 0.2 0.4 0.6 0.8 1 

redshift 


FIG. 16: The values of w(z) for the mock data set. The 
assumed value of w = — 1 is recovered accurately, with only a 
modest bias at the high redshift end. 
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FIG. 17: The potential energy of the dark energy V{z) in 
units of the critical density today. The zero redshift value is 
Vo = 0.63 ±0.05, whereas the expected value for the standard 
LCDM model is 0.7. 



redshift 


FIG. 18: The values of V{z) for the mock data set. The 
assumed value of 0.7 is recovered accurately, except for a small 
bias at the high redsihfts. 



redshift 


FIG. 19: The kinetic energy of the dark energy K(z) in 
units of the critical density today. The zero redshift value 
is Ko = 0.02 ± 0.03, whereas for the standard LCDM model 
the expected value is 0. 



redshift 


FIG. 20: The values of K(z) for the mock data set. Again, the 
correct value for the assumed cosmology, A = 0, is recovered 
accurately. 
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^Om = 0.3. The potential energy V{z) is flat, as ex¬ 
pected if the dark energy is a cosmological constant, and 
K{z) is very close to zero, also as expected in a LCDM 
model. In fact, all of the results obtained here are consis¬ 
tent with expectations in a standard LCDM model since 
P, /, w, K, and V all remain essentially constant for the 
redshift range for which the data allow a determination 
of each quantity. The two modified gravity models shown 
in Fig. and Fig. 0 do not describe the data as well 
as the standard LCDM model. 
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